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A oat  of  subroutinaa  haa  baan  davaloped  which  provida  completa 
definition  of  the  blast  environment  resulting  from  the  frte  air 
datonation  of  a one  kiloton  device  in  a ssa-laval  atmosphere. 
The  subroutines  provide  the  pressure,  density,  and  velocity  as 
a function  of  space  and  time  (from  1 ms  to  several  minutes). 

The  analytic  fits  aro  compared  with  results  of  hydrodynamic 
calculations  and  with  experimental  data.  Blast  parameters  are 
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1.  INTRODUCTION 


A requirement  for  the  description  of  the  blast  environment 
from  a 1-kiloton,  free-air  detonation  in  a sea-level  atmosphere 
has  been  acknowledged  for  many  years.  This  requirement  has  been 
fulfilled  by  a variety  of  diverse  calculations,  curves  and  ex- 
perimental data  (Reference  1) , Zn  each  case,  only  a partial 
answer  could  be  given,  such  as  peak  overpressure  versus  distance 
or  time  of  arrival  versus  distancs. 

One-dimensional  calculations  based  on  first  principles 
seemed  to  provide  a solution.  Probably  the  most  famous  of  these 
was  the  IBM  Problem  "M"  of  the  late  1940s  (Reference  2) . Con- 
sideraole  advances  in  the  state-of-the-art  of  hydrodynamic  cal- 
culations have  been  made  in  the  intervening  years  as  more  physics 
was  included  in  the  codes.  Even  now,  however,  differences  occur 
in  the  results  of  such  calculations  for  a variety  of  reasons, 
ranging  from  the  basic  differencing  method  used  to  the  personal 
preference  of  the  person  running  the  code. 

A Nuclear  Blast  Standard  (1  KT)  was  published  by  the  Air 
Force  Weapons  Laboratory  (AFWL)  in  1973  (Reference  3) . Having 
withstood  the  tests  of  time  with  only  minor  modifications,  the 
proposed  standard  has  been  accepted  and  is  presented  herein. 
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2.  THE  DATA  BASE 


Tha  sat  of  subroutines  presented  in  this  report  describe 
the  airblast  fron  u 1-kiloton,  free-air  detonation  at  saa  level. 

These  subroutine  • were  developed  as  a fit  to  tha  results 
of  one-  and  two-dimensional  radiation  hydrodynamic  and  pure 
hydrodynamic  calculations  in  Lagrangai n and  Eularian  coordinates 
using  first  and  second  order  differencing  methods  as  well  as 
acoustic  theory  for  the  very  lata  times  (Reference  4). 

A sufficient  number  of  varied  calculations  have  bean  made 
so  that  effects  of  zone  size/  difference  method/  or  other  pecu~ 
liarities  of  individual  calculations  have  been  eliminated  from 
the  fits.  In  addition  to  the  first  principle  calculations,  tha 
fit3  have  been  checked  with  experimental  data  where  available 
(References  1 and  5) . Several  previously  published  curves  for 
overpressure  versus  distance  have  also  been  consulted  (Refer- 
ences 1.  6,  and  7)  . 

Extrapolation  has  been  used  on  the  calculated  peak  over- 
pressures. This  extrapolation  has  reduced  the  error  due  to 
numerical  smearing  to  less  than  1 percent  over  a factor  of  ten 
in  zone  size.  We  are  thus  able  to  give  considerable  confidence 
to  the  calculated  peak  overpressures  and  arrival  times  as  a 
function  of  distance.  These  parameters  form  the  basic  informa- 
tion for  the  fits.  The  fit  to  the  peak  overpressure  versus 
distancs  curve  was  presented  in  August  1971  at  the  DASA  Land- 
Naval  Systsms  Long  Range  Planning  meeting  and  is  published  in 
the  Proceedings  (Referencs  9) . Since  that  time  some  minor 
modifications  have  been  made*  but  the  form  of  the  equation 
remains  unchanged.  The  maximum  error  between  fit  and  extra- 
polated calculation  is  approximately  5 percent  with  the  average 
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error  lass  than  3 percent  over  the  range  from  ten  metere  to 
several  kilometers.  The  maximum  error  occurs  in  the  early 
time  region  where  secondary  shocks  cause  perturbations  at  the 
shock  front. 

The  radius  versus  time  curve  is  accurate  to  1 percent 
over  the  range  from  ten  meters  to  several  kilometers  with  the 
maximum  error  (1  fercent)  occurring  between  IS  and  40  milli- 
seconds (160  and  270  meters). 

We  should  point  out  that  this  set  of  fits  ignores  secondary 
shocks  and  is  designed  to  give  accurate  answers  in  the  rela- 
tively ideal  case  of  single  shocks.  The  fits  are  designed  to 
account  for  the  added  impulse  due  to  secondary  shocks  without 
giving  a second  peak.  The  calculations  show  that  no  consequen- 
tial secondary  shocks  occur  beyond  a time  of  0.01  second. 

Thu  user  is  warned  that  the  use  of  any  "standard"  definition 
of  a blast  wave  for  pressure  greater  than  about  100  MPa  can 
load  to  errors.  For  radii  less  than  about  IS  meters  (~100  MPa) 
the  specific  weapon  characteristics  may  dominate  the  blast 
environment. 
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3.  OBJECTIVES 


The  1-KT  Standard  was  assigned  to  accomplish  ths  following 
objactivss  t 

1.  Determine  a functional  relation  that  describes  the 
waveforms  as  a function  of  radius  for  each  of  the  variables, 
overpressure,  overdensity,  and  material  velocity. 

2.  Provide  functional  relations  ror  peak  values  for  each 
of  these  variables  as  a function  of  radius. 

3.  Provide  a method  for  obtaining  the  time  of  arrival  of 
the  shook  front. 

4.  Determine  functional  relations  for  the  positive  phase 
length  as  a function  of  time. 

5.  whanevar  possible  then*  relations  are  required  to  be 
physically  meaningful  and  to  reproduce  the  predicted  asymptotic 
behavior  of  the  variables. 
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4.  WAVEFORM  DESCRIPTION 


The  following  ia  a general  statement  which  tracts  tha 
•volution  of  the  wavaform  from  strong  shock  to  a nearly  acoustic 
wava. 


4.1  Ovarprassura  wavaform 

Tha  first  objaotlva  was  to  find  a functional  relation 
that  describes  tha  ovarprassura  waveform,  i.a.,  a relation  between 
radius  and  overpressure  at  any  given  time.  This  relation  must 
necessarily  describe  the  following  overall  time  history  of  the 
waveform  development. 

a.  At  early  times  the  entire  waveform  consists  of 
positive  overpreeuuro  values  (positive  phase)  only.  The  over- 
pressure values  that  describe  the  waveform  decay  monotonically 
from  the  peak  overpressure  value  at  the  shock  radius  to  a smaller, 
but  positive,  overpressure  value  at  zero  radius.  As  time  pro- 
gresses! the  waveform  peak  overpressure  value  decreases  as  does 
the  positive  overpressure  value  at  zero  radius. 

b.  At  a time  t2  («  0.13  second)  the  waveform  decays 
in  such  a manner  that  it  assumes  a zero  overpressure  value  at 
zero  radius. 


o.  For  time  greater  than  t2  the  waveform  is  described 
by  both  negative  (negative  phase)  arid  positive  (positive  phase) 
overpressure  values.  Por  a brief  period  of  time  the  waveform 
decays  monotonically  from  its  positive  peak  overprossure  value 
to  a zero  overpressure  value  at  some  radius  (R  ) and  continues 
monotonically  to  a negative  overpressure  value  at  zero  radius. 

d.  Prom  some  time  (tl)  on.  the  strictly  monotonic 
decay  ceases  and  the  waveform  begins  to  assume  a waveform  which 
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will  be  referred  to  as  a well  defined  negative  phase  waveform. 
These  waveforms  can  be  described  as  follows*  The  waveform 
decays  monotonically  from  a peak  overpressure  (positive)  value 
through  a zero  overpressure  value  at  some  radius  R2,  and 
continues  to  a minimum  negative  overpressure  value  at  some 
radius,  RMIN/  for  radii  lass  than  and  continuing  to  zero 

radius,  the  overpressuro  remains  negative  but  monotonically 
returns  to  zero  overpressure  value  as  the  radius  approaches  zero. 

The  region  defined  in  paragraph  d abovo  may  be  character- 
ized by  five  basic  parameters* 

1.  The  peak  overpressure  (OP  ) . 

2.  The  radius  of  the  peak  (R  ) . 

kr 

3.  The  radius  of  the  point  at  which  the  overpressure  is 

2aro  (R_). 
z 

4.  The  minimum  overpressure  (OP^j) . 

3.  The  radius  of  the  minimum  overpressure  (Rj^) . 

Using  essentially  trial-and-error  methods,  we  found  :hat 
the  initial  decreasing  region  immediately  behind  the  shock 
could  be  fit  very  well  by  a hyperbola  of  the  form. 


OP  (r) 


-Jkri 

A (Rp-r) +b 


+ OP 

P 


(1) 


where  A and  B are  functions  of  time.  This  form  holds  from 
r ■ Rp  through  r ■ Rz»  It  is  then  necessary  to  modify  the 
function  so  it  passes  through  the  points  and  and 

returns  to  zero.  This  is  accomplished  by  multiplying  the 
hyperbola  by  the  S shaped  function  (Reference  9) 


i 

! 


,3 


J 


k 


G (r)  - 1-b 


(2) 


where  b,  c,  and  n art  functions  of  tha  basic  parameters  above. 

Unfortunately  G(r)  approaches  1 u r ♦ • but  may  be  significantly 

different  from  1 at  r ■ R . G was  modified  to  the  form 

P 


i R_-r  \ / r-R„ 

"(t> * 0(rl  (i^J  + (57*- 


(3) 


for  r.  1 r s R.  and  H(r)  ■ G(r)  for  r < R_. 
z p • 

The  expression  for  overpressure  versus  radius  at  any  time 
greater  than  0.9S  second  for  Rz  £ r ? Rp  is  given  by 


op(r>  ' U'i'v^s  * % 


and  for  r i R,  is  given  by 


^ rn\/R"r\  r-R  1 
l-bc  -2—  + ^-4- 

Av*V 


p * 


(4) 


op  (r) 


R -r 

A_(Rp-r)+8  + 0PP I 


1-b 


(5) 


In  the  time  region  described  earlier  in  this  secioh  (t  < 

0.1  second)  there  is  no  point  at  which  the  overpressure  is  zero. 
R^^  is  then  taken  to  bo  zero  and  OPMIN  is  the  overpressure  at 
(or  near)  Rj^.  The  overpressure  waveform  for  t < 0.1  second 
is  of  the  form 


OP (r)  - Becr  + OP 


MIN 


(6) 
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where  c is  a function  of  time,  and  B is  determined  by  the  peak 
conditions.  In  the  time  interval  0.1  < t < 0.95  a smooth 
switching  function  is  used  to  transform  the  early-time  function 
to  that  of  lata  time.  The  overpressure  waveform  for  any  time 
was  now  defined  in  terms  of  the  five  basic  parameters. 

The  peak  overpressure  versus  radius  curve  is  a slightly 
modified  version  of  the  formula  presented  at  the  DMA  Long-Range 
Planning  Meeting  in  1971  (Reference  8) . 

The  formula  demonstrates  our  objective  of  finding  physi- 
cally meaningful  relations.  The  formula  is  of  the  form 


OP  (R)  ■ ~ + 

p R1  R1 


R in* 


3 exp 


3 \MJ 


where  A,  B,  C,  and  R0  are  constants,  and  R is  the  shock  radius. 

The  first  term  reflects  the  reduction  in  pressure  at  early 
times  due  almost  entirely  to  the  increasing  volume.  We  realize 
that  radiation  losses  reduce  the  pressure  oven  more  rapidly 
than  1/R]  at  very  early  times,  however,  comparison  with  radia- 
tion-hydrocodes substantiates  the  inverse  R*  behavior  in  the 
region  of  the  fireball  even  prior  to  hydrodynamic  shock  formation. 

The  second  term  is  the  "normal"  spherical  divergence  term, 
the  pressure  falling  inversely  as  the  surface  area. 

The  third  term  is  a modification  of  the  asymptotic  form, 
the  shock  wave  decaying  toward  the  behavior  of  a sound  wave. 

The  asymptotic  form  before  modification  is 


OPp(R)  . — 


This  form  is  not  defined  for  R < RO  and  rather  than  using 
a switching  function,  the  term  was  modified  (see  Equation  (7)) 
to  retain  definition  for  all  R. 

The  radius  versus  time  curve  is  a combination  of  two 
regions.  At  times  less  than  0.21  second,  a function  of  the 
form 

R(t)  - at0,371  £l+ (bt+c)  • ^l-e“dt  ^ 

where  a,  b,  c,  and  d are  constants. 

For  times  greater  than  0.21  second,  an  iterative  procedure 
was  used  in  the  original  1-KT  Standard. 

The  zero-crossing  point  (Rz)  and  the  positive-phase  dura- 
tion (R+)  have  well  defined  asymptotic  forms.  These  are  given 
by 

R2  * CQt  + a (10) 

and 

R+  * b ('"  to)*  IU> 

where  a and  b are  constants,  C is  ambient  sound  soeed,  and 
RO  is  the  same  as  in  the  overpressure  versus  distance  expression. 
The  same  modification  to  the  R+  relation  must  be  made  to  ensure 
definition  for  all  R.  The  peak  radius  is  defined  as  \ 
thus  retaining  the  asymptotic  form  for  Rp.  However,  R+  is 
defined  in  terms  of  Rz,  thus,  the  necessity  of  the  iteration. 

The  iterative  procedure  was  time  consuming  and  comparison  with 
early-time  hydrodynamic  calculations  indicated  that  the  asymptotic 
form  was  not  really  applicable  to  the  early  time  (0.1  to  0.4 
second)  behavior. 
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To  eliminate  tha  necessity  of  an  iteration,  a function 
was  naadad  to  dafina  tha  poaitiva  phasa  length  (R+)  aa  a 
function  of  tima.  Tha  caviaad  1-KT  Standard  included  auch  a 
function#  however  it  aufferad  from  larger  inaccuraciea  than 
did  tha  iterative  achema  it  replaced.  The  current  Standard 
recommenda  a new  function  for  times  greater  than  0.21  second. 
The  accuracy  has  been  checked  from  0.21  to  200  seconds  and 
haa  a smaller  greatest  relative  error  than  either  of  the  two 
previous  fits.  Tha  fit  is  of  the  form 

R+  - a + b An(t)  + ctd  (12) 

where  a#  b,  c,  and  d are  constants  and  t is  time. 

A smooth-switching  function  is  used  in  the  region  0.21  to 

0.28  second.  The  point  of  zero  overpressure  is  of  the  form 

R„ ( t)  m (l-btc)  (c  t+d)  for  times  > 0.95  second,  where  e ia 
z 0 o 

the  ambient  sound  speed,  t is  time,  and  b,  c,  d are  constants. 
The  radius  at  which  the  minimum  occurs  is  given  by 

RMIj|<t)  - tt2(t)  - atb  (13) 

for  t > 0.21  second,  where  a and  b are  constants.  The  over- 
pressure minimum  ia  defined  in  terms  of  tha  peak  overpressure 
and  tima. 


4.2  Velocity  Waveform 

The  general  description  of  the  evolution  of  the 
velocity  waveform  is  similar  to  that  given  for  the  overpressure 
waveform.  However,  significant  timing  and  shape  differences 
must  be  reflected  in  the  fits. 

The  five  basic  parameters  used  to  describe  the  elocity 
waveform  arei 


1.  The  peak  velocity. 

2.  Tha  radius  of  tho  paak  value! ty, 

3.  Tha  saro  valocity  point. 

4.  Tha  minimum  valocity. 

5.  Tha  radius  of  tha  minimum  velocity. 

Tha  paak  valoeity  ia  obtained  through  tho  Rankine-Hugoniot 
relations  using  a variable  y aquation  of  stats  for  air. 

The  radius  of  the  shook  front  (tho  paak  velocity)  is  tha 
same  as  that  for  tha  overpressure  and  this  value  is  used. 

Tho  aquation  for  tha  radius  of  tha  saro  orossing  is 
similar  in  form  to  that  for  tha  overpressure  and  differs  only 
in  tho  constants. 

Tha  radius  of  tha  minimum  and  tha  minimum  velocity  also 
have  forms  similar  to  those  of  tha  overpressure. 

The  radius  of  zero  valocity  becomes  defined  at  an  earlier 
time  than  that  of  tha  overpressure  (about  85  maters) . Tha 
switch  to  the  late-time  form  takes  place  at  an  earlier  time 
(0.7  second)  and  the  form  of  the  early  time  waveform  is 

V(r)  ■ Vpk  (r/Rp)3  (14) 

where  Vpk  and  Rp  are  the  peak  velocity  and  radius,  respectively, 
and  a is  a function  of  time. 

This  form,  as  in  the  overpressure  fit,  does  not  account 
for  secondary  shocks.  The  neglect  of  velocity  fits  at  early 
times  and  small  radii  (r  < Rp/2)  can  ba  justified  by  comparing 
the  momentum  from  the  fit  and  calculation.  Although  the  veloc- 
ities may  be  significant,  the  density  is  small,  and  therefore, 
the  momentum,  energy  and  mass  flux  are  also  small  (by  two  orders 
of  magnitude)  and  may  ba  ignored. 
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A smooth  switching  function  brings  tha  velocity  waveform 
from  tha  aarly  to  tha  lata  time  form  between  0.25  and  0.7  second. 


4.3  Overdensity  Waveform 

The  overdensity  waveform  differs  from  that  of  the 
overpressure  and  velocity  in  that  it  has  a sero  crossing  even 
at  very  aarly  times  (due  to  conservation  of  mass) . The  over- 
density  has  the  following  evolution! 

a.  The  Monotonic  Decreasing  Phase.  Hare  the  over- 


density drops  from  some  peak  value  at  the  shock  front  to  a 
minimum  value  (negative)  and  retrains  nearly  constant  to  the 
canter. 

b.  The  Break-Awav  Phase.  The  shook  begins  to  separate 
from  the  hot  underder.se  fireball.  Tha  overdensity  decreases 
from  the  peak,  begins  to  level  off,  and  then  rapidly  decreases 
to  a minimum  value  where  it  remains  nearly  constant  to  the 
center.  This  nearly  constant  region  becomes  wall  defined  and 

is  refarted  to  hare  as  the  density  well. 

c.  The  Late-Time  Phase.  The  shock  is  separated  from 
the  density  well.  The  overdensity  decreases  from  the  peak  to  a 
a minimum  value,  increases  to  nearly  sero  and  than  decreases 
rapidly  into  the  density  well.  Zn  one  dimension,  this  well 
persists  for  many  seconds.  Zn  two  dimensions  it  begins  to  fill 
in  and  distort  at  a time  depending  on  the  height  of  burst. 

The  peak  overdensity  is  found  from  the  peak  overpressure 
using  the  Rankine-Hugoniot  conditions.  The  ideal  Rankine- 
Hugoniot  conditions  specify  the  density  as  a function  of  shock 
pressure  and  ambient  conditions 


(Y^l)Pl  + (y-l)P( 

p,  (y+iJp, 
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where  subscripts  1 and  0 represent  the  shock  and  ambient  condi- 
tions, respectively , this  relation  is  valid  for  variable  Y 
gasee  if  an  appropriate  average  Y is  chosen,  for  a gas  with 
specific  heat,  a linear  function  of  temperature,  the  appropriat 
average  is  given  by 


(16) 


where  Y#  is  the  average  value,  YQ  is  the  ambient,  and  y^  ie  the 
value  immediately  behind  the  shook  point. 


2<y0-1>  (Yi-1)  2{Yl-l) 

+ 1 ■ — — - + 1 


(Yfl-U-MYi-l) 


1+' 


Yi"l 


(17) 


For  ambient  sea  level  .>ir  Yq  * V5  • 1.4,  then 

2(Yl-l) 

Ye  " 1+5/2 (Yl-1)  * 1 


(18) 


This  is  the  form  used  in  finding  shock  front  conditions. 
The  value  for  Yl  is  found  by  an  iterative  procedure  through  a 
real  air  equation  of  state. 

We  recognise  that  tne  combination  of  aquations  15  and  18 
la  not  the  usual  formulation  of  the  Rankine-Hugoniot  equations 
for  density,  however,  this  form  is  in  agreement  with  available 
radiation  hydrodynamic  calculations. 


If  the  usual  equation 


is  used 


(19) 


P 


p /P 
r o 


T7 


tha  shock  dansitiss  resulting  srs  as  much  as  30  percent  greater 
than  any  current  radiation  hydrocode  results  Cor  radii  less 
than  40  motors.  Beyond  40  meters  the  two  formulations  ere  in 
good  agreement. 

The  discrepancy  between  radiation  hydrodynamic  calculational 
results  and  the  ideal  Rankine-Hugoniot  conditions  above  6 MPa 
has  not  been  resolved.  It  may  be  some  physical  process  or 
processes  or  it  may  be  the  inability  of  hydrocodes  to  resolve 
high  pressure  shock  waves. 

The  radios  of  the  pea''  is  given  by  the  radius  of  the  shock 
front.  The  radius  of  th<  is.1.0  crossing  has  the  same  late  time 
and  asymptotic  form  as  the  overpressure.  At  early  times  the 
zero  crossing  has  the  form 

Ra(t>  - atb  (20) 

This  form  is  used  for  times  less  than  0.2SS  second  after 
which  time  the  late  time  for  is  ued. 

During  the  monotonic  decreasing  phase;  the  overdensity 
waveform  is  given  by 

OD(r)  ■ A 4-  Becr  (21) 

where  A.  B.  and  c are  functions  of  time. 


IS 
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Tha  late  time  fit  hat  the  tan*  form  at  th«  latt  timt 
ovarpratturt  and  velocity  fita.  Thia  maant  that  tha  long 
lasting  danaity  wall  ia  not  daacribad  for  timet  graatar  than 
0.2  aacond.  Thia  la  alao  trua  of  two-dimensional  calculationa 
whara  tha  duration  of  tha  danaity  weJl  la  a function  of  the 
height  of  burat.  Therefore,  thia  aeama  to  be  a aufficiant 
description  of  the  ovardensity  for  bleat  environment  baaed 
on  one  dimensional  behavior. 

Tha  switching  from  early  time  to  late  time  form  takas 
place  between  0.2  and  1 second. 


* 

f 


a 
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5.  SUBROUTINE  DEFINITION 


Considerable  car*  has  been  takan  to  *naur*  that  tha  routines 
will  b*  compatible  with  a large  variety  of  computers  and  compilers. 
The  routines  are  written  in  standard  FORTRAN.  Word  length  or 
exponent  sirs  limits  should  not  be  of  oonoern  to  any  user. 

SUBROUTINE  PEAX  calculates  all  the  information  needed  for 
the  three  waveform  fits— overpressure  (OP) , overdensity  (OD) , 
and  velocity  (V)--by  a continuous  series  of  transfers  to  par- 
ticular subroutines  and  functions.  The  subroutine  is  called 
with  time  and  radius,  and  returns  the  peak  radius,  the  peak  and 
minimum  overpressure,  overdensities,  and  velocities,  the  radii 
of  zero,  and  minimum  OP,  CD,  and  V.  The  calculated  values  are 
carried  from  one  routine  to  another  through  tha  labeled  common 
block,  WFRT,  or  may  be  obtained  through  an  argument  list  ( iee 
Section  6) . 

FUNCTIONS  WFZR  .(waveform  overpressure  zero  radius) , WFDZR 
(waveform  density  zero  radius) , and  WFVZR  (waveform  velocity 
zero  radius)  calculate  the  waveform  radii  at  zero  overpressure, 
overdensity,  and  velocity,  respectively.  This  is  ths  radius 
which  separates  the  negative-positive  phase  portions  of  the 
waveform  at  the  specified  time  t. 

FUNCTION  WFPR2  (waveform  peak  radius)  calculates  the 
radius  of  the  shock  front  at  the  specified  time  (t).  The  radius 
of  the  shock  front  must  be  calculated  before  the  OP,  OD,  or  v 
peaks  can  be  determined.  The  overpressure  zero  point  must  be 
calculated  before  WFPR  is  called. 

FUNCTIONS  WFPKOP  (waveform  peak  overpressure) , WFPKOD 
(waveform  peak  over  density) , and  WFPKV  (waveform  peak  velocity) 
are  callable  with  the  peak  radius  and  they  return  the  waveform 
peak  OP,  OD,  and  V,  respectively.  The  routines  must  be  executed 

13 


in  the  above  order , a*  the  OP  peak  is  needed  in  the  OD  peak 
calculation#  and  both  are  necessary  to  determine  the  V peak. 

SUBROUTINES  WPPRMT#  WFDRMT # and  WFVRMT  are  callable  with 
a specified  time  (t)  and  radius  (r)  # and  calculate  the  0? , OD# 
and  V at  r.  These  routines  require  peak  radius#  peak  and 
minimum  parameters#  and  sero  crossings  before  being  called. 

SUBROUTINE  AIR  is  the  equation  of  state  for  air  used  at 
the  Air  Force  Weapons  Laboratory  and  is  included  as  part  of 
this  package.  It  is  used  in  evaluating  the  Rankine-Hugoniot 
relatione  for  variable  V.  The  calling  sequence  is 

CALL  AI R ( E # RKO  # GMONE ) 

where  E is  the  energy  density  in  ergs/gm#  RHO  is  the  density 
in  gm/oc#  and  GMONE  is  (yl)  • 

SUBROUTINE  MATM62  is  a model  atmosphere.  Specifically# 
it  is  the  1962  U.S.  Standard  temperate  atmosphere  with  extension 
to  700  kilometers.  The  subroutine  is  called  with  an  altitude 
(in  meters)  and  it  returns  the  atmospheric  pressure#  density 
and  temperature  at  that  altitude. 

SUBROUTINE  SCALKT  calculates  the  modified  SACH's  scaling 
coefficients  using  the  model  atmosphere  MATM62 . The  routine 
is  called  with  the  height  of  the  point  of  interest  in  m and 
the  yield  of  the  burst  in  KT.  The  routine  then  calculates 
scaling  factors  for  velocity#  density#  time#  distance  and 
pressure. 

SUBROUTINE  WELL  defines  the  fireball  density  well  to  a 
time  of  1 second.  The  fireball  is  filled  between  1 and  1.2 
seconds  because#  in  the  real  case,  a fireball  would  have  risen 
from  the  point  of  detonation.  For  a more  realistic  modal  of 
the  fireball  a three-dimensional  model  must  be  used. 
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SUBRO'JTINE  SHOCK  it  called  with  an  arbitrary  yield  and 
altitude.  Scaling  ie  done  interior  to  the  routine  and  the 
overpressure,  overdeneity  and  velocity  for  the  given  yield, 
altitude,  time,  And  radius  ara  returned. 
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6.  USE  OF  THE  ROUTINES 


For  a l-kiloton,  sea-level,  free-air  burst  PEAK  will  he 
the  only  routine  called,  Given  a time  and  a radius,  subroutine 
PEAK  returns  the  shook  radius,  the  overpra/aure,  overdensity, 
and  velocity  at  the  shook  front  and  at  the  given  radius. 

The  routines  are  muoh  faster  if  several  radii  are  called 
for  a given  time  rather  than  sever*!  times  at  a given  radius. 

The  first  evaluation  at  & nev9  time  requires  approximately 
four  times  as  much  computation  ag  subsequent  radii  at  the  same 
time. 

The  aalling  sequence  is  CALL  PEAK  (Tr  R,  SR,  OPK,  ODK,  OPR, 
ODR,  VPK,  VR)  where 


T 

is 

timq  in  seoonds 

R 

is 

radius  in  m 

SR 

is 

the  shock  radius  at  T in  m 

OPK 

is 

the  overpressure  at  SR  in  Pascals 

ODK 

is 

the  overdensity  at  SR  in  kg/m3 

OPR 

is 

the  overpressure  at  R in  Pascals 

ODR 

is 

the  overdensity  at  R in  kg/m3 

VPK 

is 

the  velocity  at  SR  in  m/s 

VR 

is 

the  velocity  at  R in  m/s 

The  results  are  for  a standard  nuclear  l-kiloton  detonation  in 

a free-air,  sea-level* environment. 

* 

For  arbitrary  yields  and  altitudes,  subroutine  SHOCK  may 
be  used. 

The  calling  sequence  is  CALL  SHOCK  (YIELD,  ALT,  T,  R,  OPR, 
ODR,  VR)  where  the  input  parameters  aret 


i 
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YIELD  - the  yield  in  kilotons  j 

ALT  - the  altitude  of  the  point  of  interest 
in  meters 

T - the  time  in  seconds 

R - the  radius  from  burst  point  to  the 
point  of  interest  in  meters 

The  calculated  parameters  are: 

OPR  - overpressure  in  pascals  at  the  point 
ODR  - over density  in  kg/m3  at  the  point 
VR  - velocity  in  m/s  at  the  point 

SUBROUTINE  MATM62  provides  ambient  atmospheric  conditions 
as  a function  of  altitude.  It  is  based  on  the  1962  U.S.  Standard 
temperate  atmosphere.  The  routine  has  five  arguments,  the  first 
is  the  altitude  of  interest  (in  meters).  The  routine  returns 
in  order,  the  pressure,  sound  speed,  density  and  temperature. 

Caution  should  be  used  whan  asking  for  atmospheric  parameters 
above  100-km  altitude.  MATM62  assumes  the  constituency  of  the 
atmosphere  to  be  tha.t  of  sea  level,  and  does  not  account  for 
the  molecular  dissociation  in  the  ionosphere. 

CALL  MATM62  (ALT,WSP,CS,WSR,WST) 

W3F  - Pressure  (pascals) 

CS  - Sound  Speed  (m/s) 

WSR  - Density  (kg/m3) 

WST  - Temperature  (°K) 

SUBROUTINE  AIR  is  the  Doan  Nickel  equation  of  state  for 

air  (Reference  10).  The  units  for  this  routine  are  cga. 

There  are  three  arguments  for  SUBROUTINE  AIR.  The  specific 

energy  and  density  are  the  first  two  arguments,  and  the  value 

of  Y-l  (wh.ire  Y is  the  ratio  of  specific  heats)  is  returned 

as  the  third  argument.  This  routine  is  valid  for  densities 

from  10  times  ambient  to  lO-7  of  ambient  and  for  energies  up 
12 

to  2 x 10  ergs/gm  above  which  it  dogrades  gracefully. 


r/ 
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EEE  - Specific  Energy  (ergs/gm) 

RRR  - Density  (gm/cc) 

GMONE  - (Y-l) , where  Y is  the  ratio  of  specific 
heats 

Occasionally  it  may  be  desirable  to  determine  the  shock- 
front  pressure  at  some  distance  from  the  burst.  Function  WFPKOP 
(waveform  peak  overpressure)  may  be  used,  however,  the  units 
are  c.g.s..  WFPKOP  returns  the  shock-front  overpressure  in 
dynes  per  square  centimeter  for  a range  R in  cm. 

The  WFPR  2 function  returns  the  shock  radius  in  cm  at  a 
time  (t)  in  seconds. 

SUBROUTINE  SKALKT  calculates  the  various  scaling  coeffi- 
cients for  arbitrary  yield  and  altitude.  The  routine  is  called 
with  the  height  of  the  field  point  and  the  yield  of  the  burst. 

HFPT  - Height  of  the  field  point  (m) 

WB  - Yield  of  the  burst  (kt) 

The  output  parameters  are  the  velocity,  density,  time,  distance, 
and  pressure  scaling  factors.  CAUTION,  the  time  scaling  factor 
is  the  Inverse. 

VS C ALE  - Scales  1-kiloton,  sea-level  velocity  to 
WB,  HFPT 

DSCALE  - Scales  1-kiloton,  sea-level  density  to 
WB,  HFPT 

TSCALE  - Scales  HB,  HFPT  time  to  1-kiloton  sea  level 

C SCALE  - Scales  1-kiloton,  sea-level  distance 
to  WB,  HFPT 

PSCALE  - Scales  1-kiloton,  sea-level  pressure  to 
WB,  HFPT 

For  implementation  and  proper  use  of  the  scaling  factors,  refer 
to  subroutine  SHOCK.  It  is  not  recommended  that  other  routines 
or  functions  be  used  independently. 
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7.  RECENT  DEVELOPMENTS  AND  CONCLUSIONS 

As  work  on  this  modal  progressed,  two  deficiencies  became 
apparent.  Several  potential  users  have  expressed  the  desire 
to  have  the  density  of  the  fireball  more  well  defined  at  lata 
times.  This  becomes  a two-dimensional  problem  because  the 
fireball  rises  in  time  while  the  shock  remains  approximately 
spherical.  The  structure  within  the  fireball  becomes  even  more 
complex  if  we  include  the  effects  of  reflected  shocks  and  shock 
torusing.  Although  these  effects  are  not  included  in  this  model, 
a fully  three-dimensional  modal  (LAMB)  is  presently  available 
to  DNA  users  through  Kaman  TEMPO,  Albuquerque. 

Four  changes  have  been  made  to  the  1-KT  Standard  routines, 
of  References  3 and  11.  The  first  is  in  the  function  WFPKOP, 
which  calculates  the  peak  overpressure  as  a funotian  of  radius. 
The  constants  in  tha  data  statement  have  been  changed  back  to 
those  of  Reference  3.  The  constants  are  as  follows t 


Parameters 

From 

TO 

AC 

3.18E18 

3.04E18 

AQ 

1.00E14 

1.13E14 

ASTAR 

9.0E9 

7.9E9 

RSTAR 

4.454E4 

Unchanged 

The  major  effect  is  in  the  region  of  peak  pressures  below 
S psi  (30  KPa) . The  new  constants  yield  a lower  overpressure 
in  this  region  rha  higher  pressure  region  is  affected  very 
little. 

The  radius  of  the  shock  front  as  a function  of  time  is 
calculated  in  the  function  WFPR2 . The  changes  in  this  routine 
affect  answers  for  times  greater  than  0.1  seconds.  Two  separate 
equations  are  used  for  times  less  than  0.21  seconds  or  greater 
than  0.28  seconds.  A linear  interpolation  of  the  two  equations 
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it  used  between  thsse  two  timtt.  Agremaent  with  the  0.5-m  SAP 
calculation  it  vary  good  for  all  timas  at  shown  in  Appendix  Z. 

A third  change  wat  made  in  tha  function  routina  WELL.  A 
constant  in  tha  aquation  for  DEPTH  waa  changed  from  -1.175E-3 
to  -1.21E-3.  This  allowt  tha  overdensitv  inaida  tha  firaball 
to  fall  to  a minimum  valua  of  -1.21  kg/nr . Examination  of 
result!  of  tha  ona-dimansional  FAB  calculation  mada  by  Kaman 
AviDyna  and  of  tha  Los  Alamos  National  Laboratory  (LAND  fire- 
ball calculation  reveal  that  tha  overdenaity  may  fall  at  low 
as  -1.22  kg/m3.  A reasonable  average  waa  taken  to  be  -1.21 
kg/m3.  Comparison  with  tha  0.3  meter  ona-dimansional  SAP 
calculation  are  alto  in  good  agreement. 

The  finaJ  changa  waa  mada  to  tha  internal  units  conversion. 
Tha  routines  are  now  programmed  for  input  and  output  in  MRS 
units.  Tha  routines  internally  convert  to  cgs  units  before 
calculations  are  mada  and  then  back  to  MKS  units  before  answers 
ara  output. 
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APPENDIX  Z 


COMPARISON  OF  THE  FIT  WITH  THE  0.50  METER  SAP  CALCULATION 

This  appendix  contain*  plot*  of  hydrodynamic  varitbl** 

« function  of  radius  at  given  timas.  Also  on  aaeh  plot  ara  tha 
rasults  of  tha  AFWL  SAP  hydrodynamic  coda  calculation  which  had 
0.5  matar  sonaa  in  tha  atmosphere.  Tha  1-KT  Standard  is  a fit 
to  many  calculations*  the  0.5  matar  SAP  calculation  is  repre- 
sentative of  those  calculations*  but  does  not  constitute  a 
"bast  fit* . Tha  SAP  data  ara  represented  by  tha  Symbol  X*  tha 
continuous  curve  is  tha  fit. 

The  radius  of  tha  shock  front  plotted  hare  was  calculated 
by  tha  non-itarativa  procedure  described  as  optional  in  tha 
main  text. 


i 

« 

I 

t 


VERPRES5URE 


RftDIUS  n 
STANDARD 


ERPRESSURE 


VELOCITY 


o.o  a.0  16.0  21.0  32.0  li.C  *8.0  56.0  6H.0 

RflDIUS  H 

SflP  IKT  STfiNDRRO  DX-.5  H 


0.0  e.0  16. G 21.0  32.0  it.O  16.0 

Rfiojus  n 

SAP  IKT  STANDARD  DX-.5  H 


KINETIC.  ENERGY 


39 


OVERPRESSURE 


40 


RADIUS  H 
STANDARD 


RADIUS  H 
STANDARD 


0 01 


0.0  8.0  16.0 


s 


4] 


0.0  0.0  16.0  21.0  32.0  11.0  18.0  36.0 

RADIUS  n 

Sfi P LKT  STANDARD  DX-.5  H 


IC  ENERGY 


&• 


KINETIC  ENERGY 


0.0  ll.o  20.0  31.0  1*3.0  50.0  60.0  70.0  60.0 

RADIUS  n 

SAP  IKT  STANDARD  DX-.5  M 


MOHENTUf 


RR01U5  H 
STANDARD 


MOMENTUM 


RADIUS  H 


0.0  20.0  10.0 


OVERPRESSURE 


OVERDENS IT 


0.0  ta.Q  60.0  I2C.0  160.0 

RADIUS  H 
SAP  IKT  STfWDARO 


Tire  - 2.000*10 


KINETIC  ENERGY 


RADIUS  n 


OVERPRESSURE 


MOMENTUM 


iS.O  60.0  120.0  160.0  200.0  240.0  230.9  320.0 

RflDIUS  H 

SfiP  IKT  STfiNOARD  DX-.5  H 


k I NET I 0 ENERGY 


0.0  18.6  60.0  120.0  160.0  200.0  210.0  200.0  320.0 

RHDJIG  r# 

SflP  IKf  SlfirtUMRO  DX-.5  h 


VELOCITY 


RADIUS  n 

SfiP  IKT  STANDflRO  DX-.5  H 


OVERDENS I TY 


RADIUS  n 

STANDARD  DX-,5  H 


MOMENTUM 


(I NET I C ENERGY 


RADIUS  N 


VELOCITY 


rroius  n 


OVERDENSITY 


RADIUS  H 
SAP  IKT  STANDARD 


u s -xo 


KINETIC  ENERGY 


RADIUS  H 
SAP  HOT  STANDARD 


VELOCITY 


03S/W  AH  3013 A 


il 


5f*>  IKT 


sap  IKT  SlftKOrtRO 


APPENDIX  II 

COMPARISON  OF  1-KT  STANDARD  WITH  OTHER  PITS 


This  Appendix  contain*  comparison!  of  savaral  parameters 
generated  from  the  1-KT  Standard  with  those  from  References  1, 
3,  6,  and  12.  when  parameters  were  not  available  from  these 
references*  the  1-XT  Standard  is  plotted  alone. 
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APPENDIX  ZXI 

PEAK  OVERPRESSURES  AS  A FUNCTION  OF  RADIUS 


Thia  Appendix  contain#  tha  peak  ovarpraaauras  aa  a function 
of  radiua  aa  oaiculatad  by  tha  1-KT  standard. 
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PEAK  OVERPRESSURE  AS 

A FUNCTION  OF  RADIUS 

R(M) 
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APPENDIX  IV 
FORTRAN  LISTING 


This  Appendix  oontsins  a FORTRAN  listing  of  the  sat  of 
subroutines  described  in  this  report. 
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APPENDIX  V 

PARAMETER  VERSUS  TXMP  PLOTS  AT  SELECTED  RADII 

This  Appendix  contains  parameter  vanui  tin*  plot*  a t 
selected  radii.  Tha  parameters  inoluda  overpressure*  over* 
prasaura  impulse*  velocity,  dynamic  prasaura * and  dynamic 
praaaura  impulse.  Tha  radii  wara  choaan  as  thoaa  corresponding 
t prasaura  iavais  of  10*000#  1000#  100#  10#  1#  and  0.S  psi. 
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DP  IMPULSE  VS.  TIME 
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CJP  IMPULSE  VS.  TIME 
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DYNAMIC  PRESSURE  VS.  TIME 


RADIUS  s 1.9216E+Q3  M 
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